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METHODS AND ARRANGEMENTS IN A TELECOMMUNICATIONS SYSTEM 
FIELD OF INVENTION 

The present invention relates to communication, and primarily 
to wireless communication. 

5 DESCRIPTION OF RELATED ART 

In a radio communication system, a bitstream or a symbol stream 
is sent on a radio channel from a transmitter to a receiver. 
Physical limitations of wireless channels, such as bandwith 
limitations, propagation loss, interference, multipath fading, 
10 inter symbol interference and time variance, present a 
fundamental challenge to reliable communication. Further 
challenges come from power limitation, size, and speed of 
devices that are employed within portable wireless apparatus . 

It is well known for a man skilled in the art that in most 
15 scattering environments, antenna diversity, i.e. the usage of 
multiple transmitting and/or receiving antennas, is the most 
practical technique for reducing the effect of multipath fading. 
The diversity technique involves the provision of a number of 
replicas of the transmitted signal to the receiver, some 
20 replicas being less attenuated by fading 

In a conventional serial data system, the symbols are 
transmitted sequentially on a channel, the frequency spectrum of 
each data symbol being able to occupy the entire available 
bandwith. 

25 A problem with serial transmitting methods in radio 
communication systems is that multipath reception having 
different delayed signals, due to different transmission paths, 
may result in Inter Symbol Interference, ISI . This occurs when 
signals being delayed, due to e.g. reflections, and having a 

30 non-negligible strength, are received having a time difference 
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as long as a data symbol or longer. The most common solution to 
this problem is the usage of an equalizer at the receiver. Based 
on estimates of channel attenuation, delay and phase of a 
channel, the equalizer tries to compensate the received signal 
5 for impact from the channel. In the time domain, the complex 
attenuation and delay description of a channel is called channel 
impulse response. The channel impulse response also enables the 
measurement of how the received energy is scattered over time by 
calculating a delay power spectrum. The range of time of which 

10 the delay power spectrum is essentially zero is called delay 
spread. A problem in the case of single carrier in a serial 
system, i.e. when the symbols are transmitted one by one after 
each other, is that for a given delay spread, the Inter Symbol 
Interference is provided over several symbols, requiring a more 

15 complex equalizer than in systems wherein the symbols are sent 
concurrently. High complexity in the algorithms requires more 
processing at the receiver, and thus high battery consumption, 
high clock rates, larger chip area, etc... It is therefore of 
great importance to keep the complexity low. 

20 A parallel data transmission system offers possibilities for 
alleviating many of the problems encountered with serial 
systems. A parallel system is one in which several sequential 
streams of data are transmitted simultaneously, so that at any 
instant many data elements are being transmitted. In such a 

25 system, the spectrum of an individual data element normally 
occupies only a small part of the available bandwith. This is 
often denoted multi-carrier transmission. 

Furthermore, the impact of a sudden short lasting impulse 
interference, e.g. sparks, is alleviated by the fact that the 
30 disturbance is spread over multiple channels, as the disturbance 
on each symbol is reduced. Hence, due to the dividing of an 
entire channel bandwith into many narrow sub-bands, the 
frequency response over each individual sub-band is relatively 
flat. Since each sub-channel covers only a small fraction of the 



original bandwith and hence the channel is flat, an equalization 
method is potentially simpler than in a serial system having 
ISI. A simple equalization algorithm and the implementation of 
differential encoding may make it possible to totally avoid 
equalization. 

The invention is directed towards a system using Orthogonal 
frequency division multiplexing OFDM that can be simply defined 
as a form of multi channel modulation where its channel spacing 
is carefully selected so that each sub-channel is orthogonal to 
the other sub-channels. Thus parallel channels may be called 
multi-carrier . 

In orthogonal frequency division multiplexing (OFDM) , a channel 
is divided into many narrow parallel sub-channels, thereby 
increasing the time of an individual symbol and reducing or 
eliminating the inter-symbol interference (ISI) caused by the 
multipath environments. On the other hand, since the dispersive 
property of wireless channels causes frequency selective fading, 
there is a higher error probability for those sub-channels when 
deep fades occurs. Hence, techniques such as error correction 
codes and diversity have to be used to compensate for the 
frequency selectivity. 

In order to provide good channel quality estimation, a pilot 
signal, i.e. a training sequence, is used. Here follows a more 
detailed disclosure on how a possible OFDM signal is transmitted 
and corrected for channel impairments . 

In this disclosure, capital letters are used to denote frequency 
domain signals whereas lower case letters denote time domain 
signals. Furthermore, the indexes k and n are used for frequency 
and time domain indices, respectively. 

First, for the purpose of coherent detection, a pilot signal is 
transmitted. A pilot is a sequence of symbols that are defined 
and well known in both transmitter and receiver. The pilot is 



denoted P(k) and is a discrete sequence of complex symbols, 
which are indexed in frequency with the parameter k and ranges 
from 0 to N-l. Prior to the transmission, the pilot P(k) is 
transposed into the time domain by an Inverse Discrete Fourier 
Transform IDFT resulting in a complex valued discrete pilot 
sequence in the time domain denoted p(n), where n is a time 
index ranging from 0 to N-l. It should be noted that 
practically, an (Inverse) Fast Fourier Transform is used instead 
of an (Inverse) Discrete Fourier Transform, because it has lower 
algorithmic complexity and therefore requires less complex 
hardware. The pilot sequence is sent over the wireless channel 
modeled by the discrete and complex valued channel impulse 
response h(n), which performs a linear filtering function of the 
sequence p(n). The resulting received sequence is r p (n), which 
is discrete and complex valued: 

r p (n)=p(n) *h(n) +v(n), where 

v(n) is the noise added in the receiver or other ambient noise 
at the antenna input. The filtering process is identified by the 
sign * between the sequence p(n) and the discrete and complex 
channel impulse response h(n) . It is a so-called convolution. 

In an OFDM system, the signals are mainly processed in the 
frequency domain. Therefore, the received N sample sequence is 
Discrete Fourier Transformed (the opposite of Inverse Discrete 
Fourier Transform) resulting in 

R p (k)=P(k) -H(k)+V(k), where 

the sequence P(k) is a known sequence, whereas V(k) is a random 
sequence and therefore best modeled as a stochastic process. 
H(k) is a channel transfer function which is unknown prior to 
the transmission. V(k) is noise and therefore random. 

It should be noted that the convolution * is transformed into a 
multiplication under the Discrete Fourier Transformation. Due to 
implementation reasons (i.e. complexity) a Discrete Fourier 
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Transform is seldom performed but instead an equivalent lower 
complexity, so called Fast Fourier Transform, is carried out. 
The equivalent is done for Inverse Discrete Fourier Transform, 
i.e. the Inverse Fast Fourier Transform is used. 

An estimate of the function H(k), which is required for an 
equalization step, is searched. Based on the received signal 
R p (k), and prior knowledge of the sequence P(k) the function 
H(k) can be determined according to: 



R p (k) 
P{k) 



10 It should however be noted that this is only an estimate, since 
noise is present. If the noise level is low in comparison to the 
received signal of the sequence R p (k), the estimation of the 
channel is precise. On the other hand, if the noise level is 
high, the estimation is poor. 

15 When the channel has been estimated, data can be transmitted. 
The data sequence (discrete and complex valued) is denoted S (k) , 
where k is the frequency index ranging from 0 to N-l. Prior to 
the data transmission, the sequence S(k) is Inverse Discrete 
Fourier Transformed, or preferably Inverse Fast Fourier 

2 0 Transformed to a time representation s (n) , which is discrete and 
complex valued. It is assumed that the time elapsed since the 
channel estimation was performed is so short that the channel 
characteristics have not undergone any significant change. 

Accordingly, when the sequence s (n) is transmitted, it will 
25 experience the same channel characteristics as the pilot signal 
did. Of course, the noise is different at reception, which is 
here indicated with w, instead of v. 



r s (n) =s (n) *h(n) +w(n) 
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When the received signal has been Discrete Fourier Transformed, 
or preferably Fast Fourier Transformed, the following signal is 
obtained 

R s (k)= S(k) • H(k) + W(k) 
5 As the channel is known, the sequence S(k) can be estimated as 
R(k) RAk)-P{k) 



5(fe) = 



H{k) R ^ 



If desired, more data symbols can now be sent. At each time 
instant, the channel estimation is used to compensate for 
channel impairments, i.e. equalization. However, the channel 
10 characteristics are likely to change over time and therefore may 
new pilot signals be sent regularly having a fixed time 
interval . 

However, a man skilled in the art understands that the 
disclosure above is related to an extremely trivial case. In 

15 practice, this is probably not the way equalization would be 
performed. Instead, equalization is often an integral part of 
forwarding error correction decoding. Furthermore, a particular 
a priori knowledge of the channel, such as delay spread or 
statistical characteristics, can be used to improve the quality 

20 of the channel estimate which in turns improves the quality of 
the estimated data S(k). 

Transmitter diversity, i.e. the usage of at least two 
transmitter means, such as transmitting antennas, is an 
effective technique for combating fading in mobile wireless 
25 communications, especially when receiver diversity, i.e. the 
usage of at least two receiver means, such as antennas, is 
expensive or impractical . 
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The performance gain of linear transform based diversity with 
ideal maximum likelihood sequence estimation (MLSE) and an 
arbitrary number of transmitter antennas is investigated and 
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compared to receiver diversity in the article by J.H. Winters 
"The diversity gain of transmit diversity in wireless systems 
with Rayleigh fading", Proc . 1994 IEEE Int Communications Conf, 
Chicago, IL, June 1994, pages 1121-1125. 

5 The published International patent application WO 99/14871 
discloses a transmit diversity method for wireless 
communication. In one illustrative embodiment, two transmit 
antennas and a single receive antenna are used, providing the 
same diversity gain as maximum- ratio receiver combining scheme 
10 with one transmit antenna and two receive antennas . 

More recently, space time coding has been developed for high 
data rate wireless communication. Space time coding schemes are 
based on coding signals in time as well as space, i.e. over 
multiple transmit antennas. Typically, a number M of transmit 

15 antennas and a number N of receive antennas are employed. In the 
article by G. Foschini, "Layered space- time architecture for 
wireless communication in a fading environment when using multi- 
element antennas", Bell Labs Technical Journal Autumn 1996, the 
author describes a multiple transmit and receive antenna system 

20 for which the channel capacity scales linearly with the minimum 
number of employed antenna at the transmit and receive side. It 
should be noted that transmitter diversity is sometimes 
classified under the field of space time coding although merely 
a single antenna is employed at the receiving side. In the 

25 latter context, it is clear that space time coding may 
alternatively to high channel capacity instead be exploited for 
increased robustness performance under e.g. fading conditions as 
transmit diversity methods mainly address. 

In the article by D. Agarwal et . al . "Space-time coded OFDM for 
30 high data rate wireless communication over wideband channels", 
Proc. 48 th IEEE Vehicular Technology Conf., Ottawa, Canada, May 
1998, pp. 2232-2236, space- time coding with OFDM has been 
studied. However, decoding space time codes requires channel 
state information, which is usually difficult to obtain, 
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especially for time-variant channels having dispersive fading. 
Said paper assumes ideal channel state information. 

Another article, also studying space time coding in OFDM having 
assumed known channel state is A. Van Zelst et al . "Space 
division multiplexing (SDM) for OFDM systems" presented at VTC 
2000. 

in the article "On channel estimation in OFDMA systems" by Jan- 
Jaap Van de Beek et al. in proceedings IEEE Vehicular Technology 
Conf. Vol. 2, July 1995, page 815-819, methods are disclosed for 
channel estimation in OFDM with one transmitter antenna. Minimum 
Mean Square Error (MMSE) and Least Square (LS) estimators with 
reduced complexity are developed that exploit the assumption of 
finite length impulse response. They also address the issue of 
spectral leakage for a discrete time channel impulse response in 
a channel having continuous path delays. Considering the aspect 
of single transmit antenna channel estimation, the drawback of 
the proposed methods are that the algorithmic complexity could 
be reduced even further with preserved performance. 

For OFDM systems using space-time coding, two or more different 
signals are transmitted from at least two different antennas 
simultaneously. The received signal is the superposition of 
these signals, usually having the same average power. If the 
channel parameters corresponding to each transmitter and 
receiver antenna pair, in a system comprising at least two 
transmitting antenna and at least one receiving antennas are 
estimated by the approach developed previously in the article by 
Y. Li et al. "Robust channel estimation for OFDM systems with 
rapid dispersive fading channels", IEEE Trans. Commun, Vol. 46, 
pp 902-915, July 1998, the signals from other transmitter 
antenna (s) will create interference. The signal-to-interference 
ratio will always be very poor, and the MSE (Minimum Square 
Error) of the estimation will therefore be very large. Hence, 
novel-parameter estimation approaches are desired for 
transmitter diversity using space- time coding. 
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A simple channel estimation for multiple transmit antennas is 
presented in the article by N. Seshadri et al . "Two signaling 
schemes for improving the error performance of frequency 
selective division duplex transmission systems using transmitter 
5 antenna diversity" in Int. Journal of wireless information 
networks, Vol. 1, No. 1, 1994. The authors propose a method for 
channel learning where one pilot sequence is multiplexed in time 
and over at least two antennas. Hence the pilot is sent on 
antenna 1 first, antenna 2 second and so on. The shortage of 
10 this method is that resources are occupied of sending training 
sequences instead of data . 

In the article by Ye (Geoffrey) Li et al., "Channel estimation 
for OFDM systems with transmitter diversity in mobile wireless 
channels", IEEE Journal on selected areas in communications, 

15 Vol. 17, No. 3, March 1999, the parameter estimation for OFDM 
systems having transmitter diversity is disclosed in detail. In 
said paper transmitter diversity using space-time coding for 
OFDM systems is studied. The channel parameter estimation 
approach is developed. The channel parameter estimation is 

20 crucial for the decoding of space-time codes, and the MSE bound 
for these estimation approaches is derived. Therefore, OFDM 
systems having transmitter diversity using space-time coding can 
be used for highly efficient data transmission over mobile 
wireless channels. However, the deficiency of the method 

2 5 disclosed in the above article is that the proposed estimation 
technique is very complex. The structure of the transmitter and 
associated signals is not discussed at all. 

SUMMARY OF THE INVENTION 

A problem with some existing methods for channel estimation in 
30 space time coded systems is that the existing methods consume a 
lot of bandwith since a training symbol is sent for each 
t r an smi t an t enna . 
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A further problem with at least one state of the art method 
using concurrent channel training sequences transmissions is 
that a preferred estimation technique is very complex both from 
algorithmic and hardware implementation point of view. 

5 The present invention relates to a method for providing 
multiple channel estimation in the field of space time coding 
in a radio communication system. Space time coding includes 
transmitter diversity, space (-time) multiplexing and other 
complex use of signal coding in time and space. 

10 The inventive solution of the problems is to send pilot signals 
concurrently instead of consecutively in order to increase 
bandwith efficiency during a multitude of channel transfer 
functions in OFDM and doing so, guaranteeing non-interfering 
channel estimates at low processing cost for both transmitter 

15 and receiver side, in relation to prior art. 

A purpose of the invention is to provide a channel estimation 
technique, which is very bandwith efficient. 

A further purpose of the invention is to provide a channel 
estimation method having a low delay and being very processing 
20 efficient. 

A yet further purpose of the invention is to ensure orthogonal 
and independent channel estimate for up to N transmit antennas, 
the number N being the number of sample values used in the 
receiver and the N channels being frequency flat. 

25 A further purpose of the invention is to reduce hardware 
complexity when providing channel estimation in OFDM- systems. 

A yet further purpose of the invention is to attenuate unwanted 
noise . 

Another purpose of the invention is to be able to provide a low 
30 complexity and robust channel estimation method for a single 



transmitter antenna case in addition to the more general 
multiple antenna scenario. 

An advantage of the invention is that it is bandwith efficient, 
since merely one synchronization symbol time slot is required 
for multiple channel transfer estimations . 

An advantage of the invention is that it is latency and 
processing efficient, since the major parts of the processing 
of the independent channel estimates, per antenna, are 
performed concurrently. 

A yet further advantage of the invention is that orthogonality 
is implicitly ensured for up to N transmit antennas by using 
the coding by N times (Inverse) Fast Fourier Transforms 
provided that no channel delay spread exists. This is always 
true disregarding actual content in the training sequence P(K), 
which therefore may be optimized for any other purpose, e.g. 
low peak to average ratio . 

A further advantage of the invention is that hardware 
complexity is very low, since the major part of the processing 
relies on hardware functions inherent for a modem used in a 
OFDM system. 

A yet further advantage of the invention is that the channel 
estimates are very exact due to the stuffing with zeros where 
noise is prevalent that acts as a noise attenuation mechanism. 

The term "comprises /comprising" when used in this specification 
is taken to specify the presence of stated features, integers, 
steps or components but does not preclude the presence or 
addition of one or more other features, integers, steps, 
components or groups thereof . 

Further scope of applicability of the present invention will 
become apparent from the detailed description given 
hereinafter. However, it should be understood that the detailed 
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description and specific examples, while indicating preferred 
embodiments of the invention, are given by way of illustration 
only, since various changes and modifications within the spirit 
and scope of the invention will become apparent to those 
5 skilled in the art from this detailed description. 



BRIEF DESCRIPTION OF THE DRAWINGS 



The invention will now be described in more detail with 
reference to preferred exemplifying embodiments thereof and 
also with reference to the accompanying drawings, in which: 



10 Figure 1 is a schematic view illustrating pilot signals 

being sent separately typical for a state-of-the- 
art method . 

Figure 2a is a schematic view illustrating pilot signals being 
sent concurrently. 
15 Figure 2b is a schematic view illustrating pilot signals being 

sent concurrently from different nodes. 
Figure 3 is a schematic view illustrating the usage of cyclic 

rotation according to the invention. 
Figure 4 is a diagram illustrating a composite channel 
20 impulse response estimate provided with some added 

Gaussian noise. 

Figure 5 is a diagram illustrating further method according 

to the invention. 
Figure 6 is a schematic view illustrating the process at the 
25 receiving antenna. 



DETAILED DESCRIPTION 

In Fig. 1 is depicted a communication system, wherein training 
sequences are transmitted sequentially in an OFDM system 
comprising several transmitting antennas. The training 
3 0 sequences are transmitted in such a way that they do not 
overlap each other over time. The system comprises a 
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transmitter 110 provided with antennas xl, x2, x3 . In figure 1, 
though only three antennas xl, x2 , x3 are shown, a man skilled 
in the art understands that more than three antennas may be 
used. The system also comprises at least one receiver 130 
5 provided with at least one receiving antenna yl, y2 . Fig. 1 
illustrates a state-of-the-art approach of performing multiple 
channel estimates in OFDM. Several training symbols pi, p2, p3 
are sent one by one after each other from the different 
antennas xl, x2, x3 to a receiving antenna yl on channels 150, 
*lS 10 151, 152. 

t|J In Fig. 2a is illustrated a system according to the invention 

Zj as well as one state of the art method. The system provides a 

\M multitude of multiple training sequences that are transmitted 

at the same time. The OFDM system in Fig. 2a comprises a 

r* 15 transmitter 210 provided with antennas xl, x2, x3 and at least 

u 

hi one receiver 230, e.g. a wireless terminal, provided with at 

B least one antenna yl . Training sequences pi, p2 , p3 are sent 

11 concurrently from the antennas xl , x2 , x3 to the receiving 

antenna yl on channels 250, 251, 252. In contrast to the state 
20 of the art method, the invention uses training sequences pi, 

p2 , p3 having a mathematical relationship with each other 

enabling distinguishing estimates in an optimal and simple 

manner . 

In Fig. 2b is illustrated a system according to the invention. 
25 The system provides a multitude of multiple training sequences 

that are transmitted at the same time. The OFDM system in Fig. 

2b comprises a transmitter 610 provided with antennas xl, x2, 

and a transmitter 611 provided with an antenna x3, and at least 

one receiver 230, e.g. a wireless terminal, provided with at 
30 least one antenna yl, y2 . Training sequences pi, p2 , p3 are 

sent concurrently from the antennas xl , x2 , x3 to the receiving 

antenna yl on channels 250, 251, 252. 

Fig. 3 illustrates transmission of pilot signals according to 
the invention. A known channel estimate sequence P(k) in the 
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frequency domain is provided to a block 308. In the block 308, 
the sequence P(k) is Inverse Fast Fourier Transformed to a 
sequence p(n) . The sequence p(n) is fed to a block 360, to a 
block 381 and to a block 382. In the block 360, a cyclic prefix 
5 CP is inserted to precede the sequence p(n). In a further 
embodiment, a cyclic suffix could be used. The cyclic prefix CP 
mitigates intersymbol interference (ISI) effects. Then the 
sequence p(n) is provided to a Digital-to-Analog (D/A) 
converter 370 where it is converted to an analog signal. Then 

10 the D/A converted sequence p(n) is provided to a first antenna 
xl. In the block 381, the sequence p(n) is cyclically rotated 
by a predetermined step comprising a predetermined number n' of 
positions in the sequence and thereby transformed to a sequence 
p(n-n'). Thereby, the sequence p(n-n') is provided to a block 

15 3 61. In block 3 61, a cyclic prefix CP is inserted to precede 
the sequence p(n-n'). Thereafter, the sequence p(n-n') is 
provided to a Digital-to-Analog converter 371 where it is 
converted from a digital to an analog shape. Then the D/A 
converted sequence p(n-n') is provided to a second antenna x2 . 

20 The two sequences are then transmitted concurrently from the 
antennas xl and x2 to a wireless terminal 330 provided with at 
least one antenna yl, y2 . In an embodiment comprising three 
antennas, the training sequence p(n) is in the block 382 
cyclically rotated by a second predetermined step comprising 

25 (n''-n') positions in the sequence and thus transformed to a 
sequence p(n''). The sequence p(n-n'') is thus provided to a 
block 362. In block 3 62, a cyclic prefix CP is inserted to 
precede the sequence p(n-n''). Thereafter, the cyclically 
extended sequence p(n-n'') is provided to a Digital-to-Analog 

30 converter (D/A) 372 where it is converted from a digital to an 
analogue shape. Then the D/A converted sequence p(n-n'') is 
provided to a third antenna x3 . The three sequences are then 
transmitted concurrently from the at least one antenna xl, x2 
and x3 to the wireless terminal 330 provided with the antenna 

35 yl . A man skilled in the art understands that the inventive 
system may comprise more than three antennas. 
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The received discrete time signal for two transmit antennas and 
one receive antenna is: 

ri (n) =p (n) *hn (n) +p (n-n' ) *h 2 i (n) +vi (n) , where 
hn and h.21 are channel impulse responses, and 
5 Vi is a noise source associated with the receive antenna. 

A result similar to a cyclic convolution is obtained when an N 

p point portion is cut out from an OFDM symbol, i.e. the removal 

tti 

t Q of the cyclic prefix, after the last signal representing the 

jfi previous symbol has arrived. The integer N defines the number 
£j 10 of positions in a Fourier Transformation. 

J3 Then, in the discrete frequency domain, the sequence ri(n) is 

j 5S ;. by means of a Fast Fourier Transform transformed to 

k* 

Jjf R!(k)=P(k) • Hn(k)+P(k) • e <i *" n '* )/ " • H 21 (k)+ V!(k) 

1=* The sequence Ri(k) can now be divided by the known training 

15 sequence P(k) . The cyclic delay n' dependent phase term remains 
and can mathematically be associated with H2i(k) . A subsequent 
Inverse Fast Fourier transform returns a composite channel 
impulse response estimate, 

h C oitvp(n)=hii (n)+h 2 i(n-n' )+Vi(n) 

2 0 Now, a commonly used constraint in OFDM is that the cyclic 
prefix length should be slightly longer than the longest 
duration of a channel impulse, i.e. delay spread. In addition, 
the OFDM symbol duration is often chosen in order to enable the 
cyclic prefix to add maximally 20-25% overhead, but preferably 

25 less than that as both energy and efficiency is consumed. 
Provided that the number of cyclic positions shift n' between 
the two antenna paths is greater than the length of the cyclic 
prefix CP, the channel impulse responses hn and h 2 i are 
guaranteed to be extracted individually from a composite 

30 channel impulse response estimate h CO np. An example of the 
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magnitude (or real or imaginary part) in dB of a composite 
channel impulse response estimate |h CO mpl having some added 
Gaussian noise is indicated in Fig. 4. In the example 
illustrated in Fig. 4 hu and h.21 are 4 and 3 point long 
5 respectively, and N-32. 

The individual channel impulse responses are easily extracted 
from the received signal since the cyclic shift of n' positions 
is known, and the extraction can be made as illustrated by 
p examples shown in Figs. 4 and 5. The cyclic shift of n' 

•J; 10 positions for the individual channel impulse responses must be 
l|3 removed prior they are used for equalization of the channels. 

t;J Then, the channel impulse response corresponding to h 2 i is 

LB cyclically shifted n' positions in the opposite direction 

compared to the direction in which the associated pilot 
K 15 sequence originally was shifted, as this unfolds the phase 
jjjj wrapping on the corresponding channel transfer function, else 

>B being incurred by the cyclic shift. A cyclic shift is performed 

£J in a vector having N positions, the last position, or the two 

last positions being put at the beginning of the vector, and 
20 the other positions being shifted to the right, clockwise. In a 
further embodiment the positions are shifted counterclockwise. 
For each individual channel impulse response, positions being 
considered irrelevant or having non-significant energy are 
replaced with zero values, alternatively, each individual 
25 channel impulse response being more smoothly attenuated by a 
smother windowing function, see Figs. 4 and 5. A man skilled in 
the art understands that the replacement of only noise 
dependent parts of the channel impulse responses with zero 
results in a significant reduction of noise and therefore 
30 provides for excellent estimation. For each channel estimate, 
noise may be reduced with roughly 10-logio (N/delay spread) dB. 
The delay spread indicates the extent of the energy dispersion 
from different multipath components, i.e. the greater the delay 
spread, the later arrives the last multipath component. 
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Two methods for separating two channel impulses responses will 
now be described with reference to Figs. 4 and 5. The channel 
impulse response has about the same length as the cyclic 
prefix. 

5 The first method is illustrated by the diagram in fig. 4. 
According to the first method, a plurality of fixed ranges is 
defined in the discrete time domain. The number of fixed ranges 
is equal to the number of transmit antennas and are upper 
limited to floor (N/CP duration) . Floor is a mathematical 

10 function choosing the closest lower integer value of the 
argument. From each of said ranges, one channel impulse 
response is extracted. For example, in range 1, a number of 
positions have amplitudes in the composite channel impulse 
sequence significantly stronger than other positions. The 

15 strong channel impulse responses are of interest for the 
inventive method. 

As the cyclic shifts are well known to both the transmitter and 
the receiver, the receiver knows when to expect the start of 
the ranges. As the receiver also knows the length of the cyclic 

20 prefix, it also knows the end point of each range. The span of 
the range is at least the length of the cyclic prefix having 
additional extension to protect from leakage. It is assumed 
that the timing synchronization for the Fast Fourier 
Transformation window start has been performed at the preceding 

25 timing synchronization step. However, small errors in timing 
are acceptable, due to the fact that said errors are transposed 
to a small time shifts in the composite channel impulse 
response and thus each channel impulse response is kept within 
each individual range. Very large timing synchronization errors 

3 0 may require larger ranges in order to accommodate protection 
for such imperfections in preceding steps. 

It should be noted that the leakage and significant timing 
errors may rotate into the end of the composite channel impulse 
response due to cyclic properties of the Fast Fourier 
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Transformation. The leakage explains why range 1 extends 
cyclically into the end of the composite channel impulse 
response. When creating the channel impulse response 
corresponding to range 1 in the composite channel impulse 
5 response, positions of range 1 are copied, whereas the 
positions outside the range 1 are set to 0. The positions 
outside range 1 have insignificant contributions to the channel 
impulse response hn and may therefore be set to zero. 

For channel impulse response within range 2, which are the 
10 results of a rotation of the training sequence p by n' 
positions, the extraction follows the same procedure as for the 
range 1. However, the resulting channel impulse response 
sequence for range 2 is cyclically shifted back n' positions in 
such a way that most of the energy is gathered at the beginning 
15 of the channel impulse response. 

For an embodiment using more than two antennas at the 
transmitter end, the same procedure is performed, but having 
different range values and subsequent cyclic rotation. Those 
values depend on how the cyclic shift is performed for the 
20 transmission of respective cyclic rotated training sequence. 

A further method is illustrated by the diagram of Fig 5. This 
method comprises an additional step. This method employs a more 
intelligent and adaptive mechanism. For each range, a first 
position having the maximum magnitude within the fixed range is 

25 determined. Thereafter, second and third positions are 
determined, being the furthermost positions within the range on 
the respective side and away from the first position also 
having a predetermined lower threshold with respect to the 
first position. The adaptively chosen subset from the fixed 

3 0 range is copied, whereas the remaining positions are set to 
zero. This step is performed for each fixed range. The figure 
shows the case for the fixed range 2 . A cyclic rotation is 
performed in the same way as the embodiment disclosed in 
connection with Fig. 4. It should be noted that if all 
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magnitudes are within the predefined magnitude difference, the 
entire range will stay unaltered. 

A man skilled in the art understands that the above disclosed 
scheme may be used for a greater number of transmit antennas. 
5 However, at most the cyclic prefix CP, floor (N/CP duration) , 
undisturbed channel transfer function, i.e. transmission over 
channels having a minimum of interference, can be guaranteed 
per receiving antenna. In addition, some guard positions must 
be left between the individual channel impulse responses due to 
10 leakage from non-whole point delays . 

In a further embodiment of the invention, additionally 
windowing is used in order to minimize the effect of Fast 
Fourier Transform leakage when path delays are intermediately 
placed between sampling points. Leakage is an inherent 

15 characteristic of the (I)DFT process when the real, i.e. 
continuous, world is sampled and results in some interference 
leakage to neighboring sample points . One such real world 
effect is e.g. that path delays may be somewhere intermediately 
placed between the sampling instances. In one embodiment of the 

20 windowing, a window function is applied in the frequency domain 
after the division by the sequence P(k), but prior to the IFFT 
to the composite channel impulse response and subsequent 
channel impulse response partitioning. After the partitioning 
each channel estimate should be compensated for the effects of 

25 the performed windowing by performing inverse windowing. 

A well known window function is the Hanning window also called 
Raised cosine window. In a preferred embodiment, the Hanning 
windowing in the frequency domain is replaced by a 
corresponding time domain filter operation on the composite 
30 channel impulse response according to h C onp_windowed (n) =h CO mp (n) 12- 
hcomp (n-1) /4-hcomp (n+1) /4, where the delay corresponds to cyclic 
shifts of h COItp (n) . This signal processing trick allows low 
complexity implementation of the frequency domain Hanning 
window operation. 
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In an embodiment of the invention comprising four transmitting 
antennas the composite channel impulse response estimate 
becomes 

H C omp(n) =hn (n-ni' ) +h 2 i <n-n 2 ' ) +h 3 i(n-n 3 ') +h 4 i(n-n 4 ') +wi (n) 

n defines the cyclic rotation, wherein ni '^3 and 
preferably ni '= 0, rii"> ni'+ CP_duration, n 3 ' > n 2 '+ CP_duration, 
n 4 '> n 3 ' + CP_duration and n 4 '< N-CP_duration, 

It should be noted that only two Fast Fourier Transformations 
are required in block 666, fig. 6, i.e. one Inverse and one 
ordinary Fast Fourier Transformation. In a further embodiment, 
the individual extracted channel impulse responses are Fast 
Fourier Transformed to their respective frequency domain 
representation H xy (k) under the assumption that block 634 in 
fig. 6 operates with channel transfer functions. 

Practically a full Fast Fourier Transformation is used. 
However, in a further embodiment, for each individual channel 
impulse response, the last Fast Fourier Transformation is of a 
reduced type, since the major content is zeros. A traditional 
(Inverse) Fast Fourier transformation can remove some 
operations, when part of the input signals are known to 
comprise zeros, as is well known for a man skilled in the art. 
Alternatively the frequency response may be generated directly 
from the DFT definition when the number on non-zero elements in 
hconp(n) are few. 

Moreover, if the receiver exploit multiple receiving antennas, 
as is indicated in Fig. 3 the inventive process can be 
performed independently for each antenna . 

Now, a receiver structure for OFDM will be disclosed with 
reference to Fig. 6. In Fig. 6, only one antenna is 
illustrated. However, a man skilled in the art understands that 
multiple antennas may be used. Signals are received at the 
antennas yl, y2... yj, where J is an integer ranging from 1 and 
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upwards limited only by practical concerns such as space, power 
consumption, cost etc. In Fig. 6 only the first antenna yl is 
shown. The signal received at the antenna yl is provided to a 
block 612, where it is amplified to an appropriate level for 
5 reception by means of an amplifying stage, often with Automatic 
Gain Control (AGC) . Then, the amplified signal is provided to a 
block 614, being an A/D converter, where the signal is 
converted from an analog shape to a digital shape and samples 
are taken of the received amplified signal at successive times. 
,f| 10 Then, the digital signal is provided to a block 616, where 

~ frequency and OFDM timing synchronization is performed. The 

-CI 

\fi timing synchronization is often based on particular training 

symbols for this specific purpose, but other methods known for 
iJJ the man skilled in the art may be used. When the timing 

f s 15 synchronization has been performed, the signal is provided to a 
j=* block 618, where the cyclic prefix is removed. A discrete and 

complex valued sequence r P ,i(n) comprising training information 
p is provided from the block 618 to a channel estimation block 

Li. 

650, which will be disclosed in more detail below. A signal 
20 r s ,i(n) comprising modulated data from multiple transmit 
antennas and affected by the channel is provided from the block 
618 to a block 620, where a Fast Fourier transformation is 
performed. The block 620 provides a signal R s ,i(k) for each 
receiving antenna yl , y2 ... yJ per OFDM symbol interval . The 
25 signal R s ,i(k) is provided to a demodulation /decoding block 634, 
where channel equalization, de- interleaving, FEC (forward error 
correction) decoding and de-scrambling is provided. 

In order to decode the data that will be sent, the channel need 
to be estimated. The channel estimation is made in the block 

30 650. The signal r P ,i(n) from block 618 is provided to a block 
660, where a Fast Fourier Transformation is performed resulting 
in a signal R P ,i(k). The signal R P ,i(k) is provided to a block 
662, where the signal R P ,i(k) is divided by a training symbol 
P(k) . The result from block 662 is provided to a block 664, 

35 where an Inverse Fast Fourier Transformation is per formed . When 
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the training symbol sequence arrives, the time domain signal is 
switched into the channel estimation stage, where a composite 
channel impulse response h'(n) is extracted. The composite 
channel impulse response h'(n) is then fed from the block 664 
5 to a block 666. In block 666, individual channel impulse 
responses for this particular reception antenna yl and all 
transmit antennas are then extracted and fed to the 
decoding/demodulation block 634 . The decoding/ demodulation 
block 634 uses the channel estimates when the signals Rs,i(k), 
10 R s ,2(k),..., Rs,j(k) are processed. The signal processed in block 
634 is then output 630 to higher layers. 

Depending of the operation with respect to time or frequency in 
the decoding/demodulation block 634, it should be noted that 
often a frequency domain representation of the channel may be 

15 more suitable than the time domain oriented channel impulse 
response. Block 634 may for example use a frequency domain 
equalizer prior Forward Error Correction decoding. 
Alternatively, block 634 may incorporate frequency domain 
channel state information in a Viterbi decoder's path metric 

20 when convolutional coding is employed. Hence, additional FFTs 
are required for this purpose. Furthermore, said operation is 
performed in parallel for potential multiple receiver antennas. 
It should also be noted that the Fast Fourier Transform 
function in block 620, 660 and additional blocks having an FFT 

25 mechanism of the same size may use the same hardware. It should 
be noted that the disclosed invention in addition to the 
channel estimate capability also enables single channel 
estimates from a single transmitter antenna. The Inverse Fast 
Fourier function may also use the same hardware with small 

30 alterations of involved multiplication factors. 



